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We have studied numerically the statistics for electronic states ( level- spacings and participation 
ratios) from disordered graphene of finite size, described by the aspect ratio W/L and various ge- 
ometries, including finite or torroidal, chiral or achiral carbon nanotubes. Quantum chaotic Wigner 
energy level-spacing distribution is found for weak disorder, even infinitesimally small disorder for 
wide and short samples (W/L » 1), while for strong disorder Anderson localization with Pois- 
son level- statistics always sets in. Although pure graphene near the Dirac point corresponds to 
integrable ballistic statistics chaotic diffusive behavior is more common for realistic samples. 

PACS numbers: 73.22.Dj,73.63.-b,73.20.Fz,81.05.Uw 



I. INTRODUCTION 



During the last few years theoretical and experimen- 
tal interest has grown immensely on monolayer graphite 
samples, known as graphene following its pioneering fab- 
rication p], 0, 0]. This remarkable two-dimensional (2D) 
honeycomb lattice structure made of carbon atoms be- 
haves rather differently from ordinary 2D metals. Its 
most interesting quantum property is a linear dispersion 
near the Fermi energy, known as the Dirac point, where 
the conduction and valence bands of graphene touch each 
other, with the low- lying excitations behaving as massless 
Dirac fermions. Anomalous integer quantum Hall effect 
for graphene has been observed [2|, [4, which depends on 
the symmetry of introduced disorder [5] and has been dis- 
cussed in terms of relativistic Dirac theory [6[. Due to its 
extraordinary properties, such as high mobility, current- 
carrying capacity and thermal conduction, graphene is 
regarded as a good candidate to replace silicon in future 
nanoelectronics Q • The possibility of gating and process- 
ing graphene sheets into multi-terminal devices is also an 
intriguing issue for related research [8[. 

The purpose of this paper is to examine the energy 
level-statistics for finite graphene sheets (quantum dots 
in the form of nanotubes) [9[, for width W and length 
L characterized by their aspect ratio W/L. We have 
also considered various directions of graphene and various 
sorts of boundary conditions (BC), in the presence of dis- 
order. Although graphene is widely regarded as a rather 
clean material it is believed that mild disorder should 
HI IH El E,Q,E3. For example, alloy type dis- 
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order from placing at random different lattice atoms or 
vacancies in the honeycomb lattice [13], or correlated dis- 
order [14]. The latter, by changing the universality class 
from orthogonal to symplectic [14], is expected to lead to 
unusual phenomena, such as weak anti-localization with 
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metallic behavior at the Dirac point. The usual short- 
ranged disorder causes localization also in graphene [TBI , 
which can limit the performance of graphene made de- 
vices. These effects of disorder naturally reflect on the 
conduction of electrons, which can be studied by usual 
scattering methods for electron waves in graphene bil- 
liards [16[. Near the Dirac point, where the density of 
states for pure graphene vanishes, the transmission be- 
comes pseudo-diffusive in the limit (W/L » 1) inde- 
pendently of BC 0. 

Recent experimental observations for size quantiza- 
tion of electron levels in graphene quantum dots [l8[ 
led to theoretical consideration of level-statistics in tight- 
binding models [l9j]. These numerical simulations for spe- 
cial geometries of the honeycomb lattice (named Dirac 
billiards, with geometrical shapes such as weakly disor- 
dered circle or triangle) are in reasonable agreement with 
the observed level-statistics in the experiment jTsj. It 
was found that weak disorder at the edges of the sample 
can also produce quantum chaos, which is independent of 
the shape of grap hene dots. Our results agree with these 
observations |l8l [l9[ since we have also obtained Pois- 
son (integrable) level-statistics in the clean limit, which 
becomes chaotic for weak disorder, independently of ge- 
ometry. Moreover, our study which was carried out for 
all strengths of disorder finds quantum chaotic behavior 
for wide and short nanotubes (W/L » 1) even in the in- 
finitesimally small disorder limit, confirming the presence 
of pseudo-diffusive behavior in this case. In the opposite 
limit of narrow long cylinders (W/L « 1) the statis- 
tics becomes ballistic ^-function type, similar to that of 
the one-dimensional (ID) chain. In other words, for re- 
alistic samples having arbitrary small disorder the out- 
standing ballistic behavior of long nanotubes is replaced 
by diffusive chaotic behavior for short and wide graphene 
samples. 

The paper is organized as follows: In Sec. II we intro- 
duce the graphene lattices studied. They are of various 
sizes and geometries, corresponding to chiral and achi- 
ral carbon nanotubes. Their tight-binding Hamiltonian 
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in the presence of disorder is introduced. In Sec. Ill 
we show that for diagonal and off-diagonal disorder (in 
the case of large W/ L even for infinitesimally small dis- 
order) Wigner chaotic level-statistics is obtained while 
for strong disorder it becomes Poisson, independently of 
the type of graphene sheet. Our results for the partic- 
ipation ratio complement these findings by showing ex- 
tended eigenstates close to the Dirac point and confirm 
the presence of Anderson localization for strong disorder. 
Finally, in Sec. IV we summarize our main conclusions. 



II. THE TIGHT BINDING HAMILTONIAN 

The tight-binding Hamiltonian of tt bonding in 
graphene is 
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where c$(cj) annihilates (creates) an electron at the sites 
of the honeycomb lattice, S{ represents the on site energy 
which in the case of diagonal disorder is a random number 
in the range [— w/2, w/2], where w denotes the disorder. 
For off-diagonal disorder Si = and the logarithm of the 
nearest neighbor hopping element 7 takes random values 
in the range [—w/2,w/2] [20[. 

In the case of periodic BC only in one direction, 
the other being much longer with hard wall BC, the 
graphene samples represent finite single wall carbon nan- 
otubes with W/L « 10. These are effectively 
one-dimensional objects, build up by wrapping sheets 
of graphene into very long cylinders, and in modern 
technology underpin the evolution to nano-electronics. 
Their electronic properties depend on the wrapping direc- 
tion which is characterized by the so-called chiral vector 
(n,m), expressed in terms of graphene primitive unit lat- 
tice vectors 7f 1, if 2, = \~cti\ — V3cic-C — Q>, where 
ac-c is the distance between two atoms and a the hon- 
eycomb lattice constant. Since various possibilities exist 
to roll up graphene, depending on the rolling direction 
one obtains armchair, zig-zag and chiral nanotubes 
The armchair nanotubes characterized by a chiral vector 
(n,n) can be obtained by rolling the graphene sheet along 
one of the three vectors joining a honeycomb lattice site 
to its nearest-neighbors. The perimeter of the (n,n) arm- 
chair nanotubes consists of n hexagons along the rolling 
direction connected by n single bonds and they are always 
metallic. The (n,0) zig-zag nanotubes are defined by a 
rotation of the graphene sheet by 90° and the n hexagons 
followed by n single bonds lie now along the longitudinal 
axis of the tube (the rolling direction is perpendicular). 
The (n,0) nanotubes are metallic only when n is a multi- 
ple of 3. In all other cases, that is rotating the graphene 
sheet in between the previous two and then rolling it up, 



is equivalent to placing the hexagons followed by single 
bonds along intermediate directions. This third type of 
carbon nanotubes is called chiral (n,m), n 7^ m and its 
band structure is metallic only when n — m is a multiple 
of 3 0. 

Our approach consists of studying the effects of disor- 
der in graphene by obtaining numerically the eigenval- 
ues and eigenvectors of the Hamiltonian H in the pres- 
ence of disorder. In our computations we have fixed the 
short-ranged site (diagonal) and bond (off-diagonal) type 
disorder for many random samples. The statistics of the 
obtained eigen-solutions are subsequently analyzed in or- 
der to address the quantum chaos issues. Since for the 
(n,n) armchair nanotube a unit cell along its long di- 
rection consists of two slices around the tube, (with 2n 
carbon atoms each making An atoms in total the total 
number of hexagons covering the unit cell is Nh ex — 2n) . 
The number of such unit cells along a finite nanotube is 
N c so that their size can be classified in terms of n and 
N c . 

In the absence of disorder (w = 0) the dispersion of 
pure graphene sheet taken in the correct orientation, that 
of armchair nanotubes along the x-axis with their rolling 
direction along the vertical ?/-axis, is [2l| 

E(k x , k v ) = ±7 ( 1 + 4cos(^|^) cos(^) 
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with quantized values of k = (k x , k y ). For a finite arm- 
chair nanotube (periodic BC along y and hard wall BC 
along x) the eigenstates are labelled by two integers, via 



k x CL 



1 Joe 
2 



ky CI 



2tt 
V3n 



(3) 



with j x = 1,...,7V C and j y = l,...,2n. The output is 
2nN c positive and exactly e qual negative eigenvalues (a 
fact of sublattice symmetry [22[, since the two types of 
atoms A and B in graphene for nearest neighbor hop- 
ping make two interconnected A and B sublatices). In 
the case of infinite nanotubes the k x a quantization of 
Eq. (3) does not arise and the band structure of Eq. 
(2) becomes one-dimensional, a function of k x only [9[. 
For toroidal armchair nanotubes, where appropriate peri- 
odic BC are imposed also along the longitudinal direction 
k x u = jq^Jxi jx — 1? ••••} A^c- 

For the more general chiral (n,m), n 7^ m nanotubes 
the quantized component of the wave vector (n,m) along 
the perimeter of the nanotube where periodic BC are 
imposed, satisfies 
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FIG. 1: The normalized level-spacing distribution function 
P(S) from 10 samples of graphene with diagonal disorder (for 
the values of disorder strength w see figure) : (a) For graphene 
structures with a very small aspect ratio W/L « 0.04 and 
4960 sites corresponding to N c = 40 unit cells of the chiral 
(7,4) nanotube. As the disorder strength increases a gradual 
crossover is seen, from a broad ^-function (in the form of a 
Gaussian around the mean) for clean nanotubes to a Poisson 
distribution for strong disorder, (b) The same as in (a) but 
for a wide and short chiral (36,180) nanotube with N c = 1 
having aspect ratio W/L « 21 (4464 sites). In this case the 
distribution is close to Wigner even for almost zero disorder. 
(c),(d) From the (7,4) and (36,120) toroidal nanotube ge- 
ometries consisting of N c = 40, 12 unit cells corresponding to 
aspect ratios W/L w 0.04, 17.3 and sites 4960, 5760 sites, re- 
spectively. In the insets the densities of states are shown with 
the characteristic dip near the Dirac point, which remains for 
weak disorder but closes for strong disorder (w > 3), leading 
to a broad density of localized states. 



having width W = a\Jn 2 + m 2 + nrn given by their 
perimeter and Nhex — 2L 2 /(a 2 d r ) hexagons in a unit 
cell of length |T| = \[?>W/d r , d r being the greatest com- 
mon divisor of 2m+n and 2n+m. For a given value of the 
quantized k y the energy of an infinite chiral nanotube is 
a function of the continuous longitudinal component k x 
only, with — ^ < k x < +-px T [- For armchair nanotubes 

(n,n) the unit cell reduces to width W = a^n, length 
|T| = a and Nhex — 2n. Since L = N c a for N c unit 
cells the considered aspect ratio is W/L = d r /(y/3N c ) 
for chiral (n,m) nanotubes, reducing to W/L = \/3n/N c 
for armchair (n,n) nanotubes. 



III. RESULTS 
Pure graphene sheets 

Numerical computations of the eigen-solutions are re- 
quired for chiral nanotubes even in the absence of dis- 
order, since their unit cell along the tube can become 
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FIG. 2: For the same parameters as in Fig.l, but for the 
participation ratios of the corresponding eigenstates versus 
energy. The participation ratios indicate the number of sites 
where each eigenstate has a significant amplitude. For a per- 
fect extended state its value should be « 5000 sites for the 
chosen sizes. In the insets the results for pure graphene w = 
are magnified, note the linear scale on the y-axis. 

arbitrarily large depending on (n,m). For pure graphene 
we have found the expected ballictic behavior manifested 
via integrable Poisson statistics for arbitrary W/L, which 
becomes ballistic ID-like in the limit W/L « 1. This 
can be simply understood near the Dirac point where the 
density of states is linear. In this regime the integrated 
density of states is J\f(E) ex E 2 . Therefore, the statistics 
of the energy levels E should be replaced by the statis- 
tics of the unfolded squared levels E 2 which is required 
in order to make the density of states constant by fixing 
the mean level-spacing. If we use the usual quadratic dis- 
persion for the squared levels E 2 oc k 2 + k 2 the statistic 
near the Dirac point becomes equivalent to the statistics 
of integrable billiards [23] , with E 2 replaced by E from 

E j x j y = ajl+jy, Jxi j y integers, (5) 

a being an irrational number fixed by the size of the 
adopted sample. The resulting statistics for zero disor- 
der is obviously Poisson [23[ , except for very long samples 
in the ID limit where the sum of Eq. (5) should depend 
on one parameter only and the corresponding level statis- 
tics becomes a trivial ^-function. We have verified this 
behavior, in agreement with [l9j . However, for infinites- 
imally small disorder as seen in Fig. 1(b) the statistics 
changes and becomes chaotic when W / L is large (Wigner 
distribution) while it remains trivial ballistic ID-like for 
small W/L shown in Fig. 1(a) (broad ^-function). In the 
almost clean limit, the ballistic behavior for small W/L 
is replaced by chaotic diffusive for larger W/L samples. 

In the insets of Fig. 2 the behavior of the correspond- 
ing eigenstates for pure graphene can be seen on a linear 
scale, via their participation ratios. The main observa- 
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tion is that for the small W/L (Fig. 2(a) inset) the ma- 
jority of ballistic states have a higher participation ratios 
when compared to those with large W/L (Fig. 2(b) in- 
set). The distribution of the participation values becomes 
sharper for the small W/L torroidal nanotubes (Fig. 2(c) 
inset) when compared to large W/L (2(d) inset). 



Diagonal and off-diagonal disorder 



The density of states and the level-spacing distribution 
function P(S) are computed for various graphene sheets 
with aspect ratios W/L, including chiral nanotubes. The 
corresponding nearest-level spacing distribution function 
is illustrated for different values of diagonal disorder in 
Fig. 1. All the unfolded energy levels in the band are 
considered for 10 runs, making the density of states con- 
stant everywhere in the band. No significant difference 
was found when levels near the Dirac point are consid- 
ered only. For small W/L (Fig. 1(a)) the level spacing 
distribution function P(S) is shown to vary from a broad 
(^-function like curve for almost zero disorder (w = 10 -7 ) 
towards the Poisson localized limit for strong disorder 
w >>. For large W/L we observe instead (Figs. 1(b), (d)) 
a crossover not from a (5-function but from a Wigner 
chaotic distribution towards the localized Poisson limit 
in the strong disorder limit (w >>). In the insets of Fig. 
1 the corresponding average density of states is shown 
with the dip near the Dirac point, which disappears when 
the disorder increases. 

In Fig. 2 we demonstrate the participation ratio which 
measures the extend of the corresponding eigenstates on 
the lattice. Near the Dirac point the values are seen to be 
higher for small W/L (Fig. 2(a)) when compared to large 
W/L (Fig. 2(b)). This observation is related to the ob- 
served ^-function and Wigner P(S) distributions in Fig. 
1, for small and large W/L, respectively. Similar results 
are obtained for torroidal nanotubes (Fig. 2(c) and (d)). 
For strong disorder our results indicate localization of 
all the states in the band with their participation ratio 
becoming small. 

Off-diagonal disorder preserves chiral symmetry since 
it connects one sublattice to the other [20[. In Fig. 3 
we show our results for the chiral (7,4) nanotube with 
off-diagonal disorder in the small W << L nanotube 
limit. Fig 3(a) shows P(S) (with the density of states 
in the inset) and fig. 3(b) the participation ratio of the 
corresponding eigenstates. Similar results are obtained 
to the diagonal disorder case, with the exception of the 
appearance of a singularity in the density of states, which 
develops at the band center (not seen in the figure due 
to scale) [20|. 
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FIG. 3: (a) The level-spacing distribution function P(S) and 
(b)the participation ratio for off-diagonal disorder of strength 
w (see figure) for the small aspect ratio W/L « 0.04 chiral 
(7,4) nanotube geometry of N c = 40 cells (4960 sites). 



IV. DISCUSSION 



It is known that the band structure of graphene is 
metallic without a gap, having a density of states which 
becomes zero at the Dirac point. However, depending on 
the geometry of the samples a small gap can open near 
the Dirac point, e.g. for the (n,m) nanotubes which be- 
come insulating when n — m is not a multiple of 3. Our 
study of the eigensolution statistics for various geome- 
tries specified by the aspect ratio W/ L in the presence of 
disorder reveals some interesting features: First, quan- 
tum chaos becomes relevant also for weakly disordered 
graphene, with pseudo-diffusive Wigner statistics even 
for almost zero disorder when W/L is large. Second, for 
strong disorder localization occurs for all states in the 
band, in agreement with finite-size scaling transfer ma- 
trix studies at the Dirac point. Third, near the Dirac 
point, due to the linearity of density of states of pure 
graphene, a natural unfolding of levels corresponds to the 
study of squares of energies E 2 which gives integrable ir- 
rational billiard statistics. However, the pseudo-diffusive 
quantum chaotic behavior of graphene for infinitesimally 
small disorder independently of BC, for large aspect ratio 
W/L should become obvious in scattering experiments 
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from nanotubes. 

Since infinitesimal disorder is inevitable also for 
graphene the question whether ballistic or diffusive be- 
havior prevails can be answered in favor of a simple ID- 
like ballistic behavior for small Wj L and pseudo-diffusive 
for large W/L. For weak disorder quantum chaos and for 
strong disorder localization occurs, in agreement with 
other 2D disordered lattices. In conclusion, our com- 
putations for graphene in the presence of disorder show 
quantum chaos or Anderson localization for weak disor- 
der (even almost zero disorder for large W/L) or strong 
disorder, respectively. This demonstration of quantum 
chaos in graphene is in agreement with recent theoretical 
and experimental studies. 
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